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Atmospheric temperature profiles can be deduced from the 
molecular density profiles assuming hydrostatic equilibrium for the 
atmosphere where the density profiles can be obtained from the Rayleigh 
backscatter lidar signal profiles. Hauchecorne and Chanin^ have used 
this technique to measure temperature above 35 km altitude where aerosol 
scattering is considered negligible compared with Rayleigh backscattering. 
The lowest altitude at which temperature can be retrieved using a single 
wavelength lidar extends to about 30 km, depending on the upper boundary 
of the stratospheric aerosol layer. 

For the case of background aerosols, Russell and Morley proposed 
retrieving temperature in the lower stratosphere using two wavelength 
lidar signals. The technique consists of retrieving the aerosol 
backscatter profile at the long wavelength, then modeling a wavelength 
conversion factor for the aerosol backscatter at the short wavelength 
(approximately 0.35 urn). The estimated aerosol backscatter profile is 
combined with the measured signal at the short, wavelength to yield 
molecular density profile which is used to retrieve temperature. Error 

analyses of temperature retrievals from spaceborne lidar have previously 

, 2,3 

made. 

In this letter, differences between retrieved temperature T from 
lidar and atmospheric model temperature T^ are calculated for the 
following problems which except for problem b below have not been 
simulated so far. 

a. Effect of a transient thin aerosol layer at any altitude on 
temperature retrieved by a single wavelength lidar. 


b. Effect of residual aerosols which are not perfectly corrected for 
using a modeled wavelength conversion factor on temperature 
retrieved by a two wavelength lidar. 

c. Errors due to splicing two backscatter signals which are obtained 
seperately. 

d. Effect of multiple scattering. 

Fig.l shows Tidar signals normalized to the Rayleigh backscatter 
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signal expected from a model atmosphere for the above four cases. 

To calculate the temperature difference AT=T-T^, lidar parameters used 
for the LITE (Lidar In-Space Technology Experiment) simulations have been 
used. However, results denoted below can be applied to other lidar 
parameters if sufficient backscatter signals are obtained. 

Fig, 2 presents vertical profiles of the temperature difference 
AT simulated for above four cases. The small temperature difference 
at about 50 km is due to the boundary condition of temperature, T^ + 40 K, 
assumed at 90 km. The results shown in Fig. 2 are summarized as follows; 

a. At 34 km a AT of -8.3 K occurs due to the transient aerosol layer 

of scattering ratio 1.04. At 33 km it is about 1.2 K and approaches 
zero at 20 km. It is easy to distinguish a transient thin aerosol 
layer or noise with 4 percent excess signal above expected from 
Rayleigh backscatter from retrieved temperature profile. 

b. Effects of an unknown residual aerosol backscatter of < 1.2 % 
affect the retrieved temperature by < -1.25 K. As the maximum 
scattering ratio at a wavelength of 355 nm is about 1.04 (4 percent 
aerosol contribution) for background stratospheric aerosols. 



c. A discontinuity of AT equal to -4.2 K occurs at 30 km for 

a 2 percent splicing error constant below 30 km. AT approaches 
to zero exponentially downwards. To keep this temperature dis- 
continuity below 2 K, it is necessary to splice two lidar backscatter 
signals with about one percent accuracy.^ LITE will utilize only 
one photomultiplier channel and thus will not use splicing in its 
data analysis. 

d. Effects of multiple scattering is assumed to be 4,4 percent at 4 km 
and to decrease upwards in proportion to atmospheric densities.^ 

The calculated temperature difference is - 5 K at 4 km and approaches 
zero exponentially upwards. It is about - 1 K at 15 km. If lidar 
parameters such as receiver field of view and lidar system altitude 
are known, the effects of multiple scattering can be corrected to 
some extent theoretically.® 

From these simulations, it is evident that temperature can be 
retrieved for + 3 K, if Rayleigh backscatter can be measured for + 1.5 
percent. This is not unreasonable for a carefully developed lidar system. 
The LITE lidar in a 240 km orbit aboard shuttle should be able to make 
+ 3 K measurements from about 10-40 km, with vertical resolutions of 
1 km and horizontal resolutions of 300 km. 
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Figure Captions. 


Excess lidar signals normalized to Rayleigh backscatter 
signals for four cases : a. transient thin aerosol layer, 

b. residual aerosols, c. splicing error, and d. multiple 
scattering. 

Temperature difference between temperature T retrieved by 
lidar and model temperature for four cases in Fig.l. 
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Signals Normalized to Rayleigh Backscatter 
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Abstract 


An error analysis of DIAL (Differential ^sorption ^idar) measurements of 
stratospheric ozone from a Space Shuttle is discussed. A transmitter system 
consisting of a KrF excimer laser pumping gas cells of H 2 or D 2 producing output 
wavelengths in the near uv is shown to be \aseful for the measurement of ozone in 


an altitude range from 15-50 km 


I. Introduction 


Measurements of ozone in the lower stratosphere have been conducted using a 

1 2 

lidar incorporating an XeCl laser with supporting meteorological data, , 
Differential absorption lidar (DIAL) measurements of ozone in the lower 
stratosphere have been conducted using frequency-doubled dyes pumped by a Nd:YAG 
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laser# DIAL measurments of ozone in the middle and upper stratosphere above 25 
)cm altitude have been made using an XeCl laser pumping a cell of methane gas 
producing the first Stokes line of methane. In addition r DIAL measurements of 
ozone above 25 3cm have been conducted using the third harmonic of a Nd:YAG laser 
plus the fundamental output of an XeCl laser 

As above, DIAL measurements of ozone appear feasible up to an altitude of 
50 3cm by a ground-based XeCl laser. However^ it takes several hours to measure 
ozone densities at 40-45 3cm even if an efficient XeCl laser is used with an 
output energy of 250 mJ per pulse and a repetition rate of 50 Hz.^ This is due 
to the strong absorption of ozone below 40 3cm. In the DIAL measurements of 
ozone from a shuttle, this strong absorption can be avoided and global data can 
be obtained in a short time with high accuracies as denoted below. The global 
data will be useful in the analysis of ozone transport in the stratosphere and 
for studying long term trends in ozone. 

This paper will calculate the errors associated with and, therefore, 
explore the feasibility of the measurement of stratospheric ozone from a 
spaceborne instrument aboard shuttle. 
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II. Methodology of Error Analysis 

The mean ozone concentration n(z) within an altitude interval of z and 
z'=z + hz can be obtained by the DIAL measurement as given in the following; 


n(z) = 


1 
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and P^(z) (i*1,2) is the total received number of photoelectrons from an 
altitude z with a small range 6z, and Aa is the differential ozone absorption 
cross section (Ao=«a 2 - a-| ) . 3 i( 2 ) is the atmospheric backscattering 

coefficient at wavelength i and altitude z. or and are Rayleigh and Mie 
extinction coefficients, respectively. 
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In general the DIAL measurement of n(z) is made using the first term in 
Eq. ( 1 ) . For this approximation, the relative uncertainty e = An/n can be 
expressed by® 

( 4 ) 

where 

, 2 P. (z) + 2P. P. (z‘) + 2P. 1/2 

< 5 > 

2 Aa n Az i*1 P, (z) P. (z') 

•1 1 


* (B + E)/n 


is background noise, c-| denotes a statistical error due to the* signal 
and background fluctuations, and £2 ^ systematic error due to the assump- 

tion of B=E=Oe^,^ As the two wavelengths can be emitted simultaneously, the 
errors due to temporal variations of B and E ara not taken into accounte 
III. Lidar System Parameters 

Figure 1 depicts the transmitter system for the proposed DIAL system. A 
KrF laser pumping gas cell of either H 2 or D 2 produces four wavelengths; the 
fundamental of KrF plus three Stokes lines of H 2 or D 2 • ® ^ The output 
energy of the KrF laser is 300 mJ with a pulse width of 10-20 ns and a pulse 
repetition frequency of 100 Hz* The energy conversion efficiencies of these 
Stokes lines S^i , 82 ^ and S 3 are assumed to be 20, 20, and 10 percent, 

respectively* In these simulations, ozone absorption cross sections from 
13 14 

Inn and Tanaka , the U*S* Standard Atmosphere , and a mid-latitude ozone pro- 
file are used, and the aerosol models employed are taken from Shettle and 




Fenn 
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IV . simulation 


The total number of signal photoelectrons P^Cz), for a vertical range of 

1 km (5za1 km) emd the output energies shown in Tedsle I, are shown in Fig. 2* 

The total shot number of laser pulses is 10,000 shots. From the figure, it is 

clear that photoelectron counters should be used for the detection system. A 

fast speed counter and photomultiplier tube (PMT) with fast rise time having a 

counting linearity up to ^U^out 50 photoelectrons per shot for a time gate inter- 
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val of 6.66 us (6z«1 km) is required. 

Figure 3 shows the error analysis for the wavelength pairs consisting of 

4 

the KrF laser and H2 cell. The horizontal resolution is 800 km (10 shots at 

100 HZ). The vertical resolution is 1 km (5z=Az=1 km) below 50 km altitude amd 

2 km ( 6zaAz>2 km) above it. The solid lines in the figure show the relative 

uncertainties for a moonless nighttime when the total number of photoelectrons 

due to background noise is below 5 in the wavelength range of 248-360 nm* The 
dashed lines show the uncertainties for a case of higher noise-occurrence- 
probabilities of 0e1 background photoelectron per shot for a wavelength region 
of 248 - 313 nm, 0*2 background photoelectron per shot for •320nra, and 0*5 back- 
ground photoelectrons per shot for *360 nm. The kinks in the line around' 50 km 
in the figure are due to the difference of vertical resolution above and below 
50 km. 

For the wavelength pairs of 248.4 - 313.0 nm and 277.0 - 313.0 nm, the 
relative uncertainties of ozone densities measured by the DIAL technique are 
within 10 percent for 37 - 48 km altitude for a moonless nighttime. In this 
altitude range, the errors due to neglecting B and E are less than 0.3 percent. 
If the output energies at these wavelengths are increased by a factor of 2, the 
relative uncertainties are decreased by a factor of /2. These wavelength pairs 
are most useful for ozone measurements in the region of 40-45 km where the 
effect of CFM's on the long term ozone concentration is considered to be large. 


The vavelength pair of 313*0 - 359.8 nm is available for ozone measurements 
in the lower altitude range of 12-27 km. In this case, the errors in n(z) due 
to neglecting the term - cij^2 large below 18 km altitude. 

The dotted line in Fig. 3 shows the relative uncertainties of ozone densities 

calculated by using £' instead of E in Eg. (6) vdiere E'«(0*25 x 

(a^l - and n is estimated by the first term in Eg. (1) plus 

factor 0.25 corresponds to a standard deviation of 
14 

atmospheric density. 

The relative uncertainties of the DIAL measurement of ozone increase in the 
altitude range 27 - 34 km. This can be overcome by using a gas cell of D 2 
instead of H2* Figure 4 shows the relative uncertainties of the DIAL measure- 
ment of ozone using the stimulated Raman scattering lines of D 2 * The solid 
lines in the figure show the relative uncertainties for a moonless night, and 
the dashed lines show them for a case of higher noise-occurrence-probabilities 
as mentioned above. As can be seen, the wavelength pair 291.8 - 319.6 nm is 
useful for the DIAL measurement of ozone in cui altitude range of 25 - 40 km. 
Figures 3 and 4 show that if gas cells of H 2 and D 2 are pumped simultaneously by 
a high powered KrF laser, accurate DIAL measurements of ozone are possible in an 
altitude range of 15 - 50 km using these combinations of wavelengths. 

A. Effect of Aerosol 

Stratospheric aerosols will have am effect on DIAL ozone retrievals at the 
altitudes ^4iere they reside, approximately from the tropopause to 10 km above 
the tropopause. The magnitude of their effect will depend of course on the 
amount present in the stratosphere which is dependent on volcanic perturba- 
tions. Effects of the aerosol model assumed in the simulations on the total 


I 


errors for the wavelength pair of 31 3.0 - 359.8 nm have been auialyzed. Figure 5 
shows three aerosol models for 313.0 nm: background, aged volcanic, and fresh 
volcanic. The results of the uncertainties in ozone density for these models 
are shown in Fig. 6 , where E' is also used to calculate ££• fresh vol- 

canic aerosol model (dashed line), the relative uncertainty is large. The mini- 
mum of the uncertainty around 19 km is due to am assumption that the aerosols 
are constant vertically in spite of the maximum of the content. In this case, 
the following procedure might decrease the uncertainty. The aerosol distribu- 
tion is first obtained from the wavelength of 359.8 nm. Then the effects of 

aerosols on the 313.0 nm wavelength are extrapolated from 359.8 nm. These cases 

1 9 

are also encountered in the DIAL measurement of ozone in the boundary layer. 
Unfortunately, both fresh volcanic aerosols and boundary layer aerosols are 
difficult to model. Another method to decrease the errors is to use an appro- 

, 3 

priate wavelength pair^with a small difference of wavelength. 

B. Effect of temperature 

The ozone absorption cross sections are temperature dependent. Figure 7 

shows the relative variations of the absorption cross sections at seven wave- 

20 ^ 

lengths in a temperature range 214-291 K. These wavelengths are near €he 

laser wavelengths used for the DIAL simulations. In the wavelengths shorter 

than 291 .4 nm, the relative temperature variations of the cross sections are 

several percent, but they are 20-30 percent in the wavelengths longer than 

2 

300 nm. For the cross section variations, an additional term £3 is included in 
Eq. (4) , \^ere 
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(7) 

For the wavelength pairs used for the simulations so far, it was approximated 

that cuid 02 /CJ ^ 0*1 • Therefore, in spite of the large temperature 

2 2 
variations of 02 , is roughly equal to or smaller than 2(Aai/ai) • To obtain 

the ozone density with a high accuracy, therefore, it will be necessary to meas- 
ure temperature with an accuracy of 5-10 K# 

Figure 8 shows the U.S, Standard Atmosphere model of temperature, Tg# 

The arrows in the figure show the lowest and highest mean monthly temperatures 

14 

obtained for any location between the Equator and pole. Atmospheric tempera- 
tures will change greatly during a shuttle orbit covering both hemispheres. 
Temperature retrievals from the non-DIAL or backscatter lidar measurements at 

359.8 nm are obtained by using the perfect gas law and hydrostatic equili- 
21 

brium. The U.S. atmosphere model is also used for the simulation. The effect 
of aerosols is neglected above 30 km. As a boundary condition of temperature, 
two cases of Tg+40 K (solid line (1) to right of center line) and Tg - 60 K 
(solid line (2) to left) at 90 km are assumed. As shown in the figure, devia- 
tions from the model (center line) are very small below 65 km altitude for the 
two boundary conditions. At 60 km, for example, this difference is less than 

0.7K. Temperature retrievals below 60 km, therefore, do not depend on the 

2 1 

boundary conditions. These measurements are accomplished with only one 
wavelength. 




Aa, 2 Aa, 2 Ao_ 2 Aa_ 2 
_^) (— ^) + {— ^) ( —) 


1 
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Figure 9 shows the accuracy in temperature measurements for various output ' 
energies per pulse using the above one wavelength technique. The boundary 
condition of temperature at 90 3on is assumed to be Tg+40 K. The vertical 
resolutions are 1 km (5z*Aza1 km) below 50 km altitude and 2 km (Sz*Az=2 km) 
above it. As the output energy of the laser increases, measurement errors of 
temperature decrease. Even for low output energies, it appears possible to 
measure vertical distributions of temperature with an accuracy of 2-5 K in the 
region of 40 km altitude. Simulations of temperature retrievals for LITE (Lidar 
In-Space Technology Experiment), for example, using the tripled wavelength of a 
NdtYAG laser (355 nm) shows similar results in that temperatures can be 
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retrieved to ±3K between 15 and 30 km during background aerosol conditions, 

V. Conclusion 

A shuttle lidar system based on a KrF excimer laser pumping a Raman cell of 
or D 2 producing the fundamental wavelength and the Stokes lines of H 2 and D 2 
are shown to provide an accuracy useful for the DIAL measurement of ozone in the 
middle and upper stratosphere. The relative uncertainties of ozone densities 
measured at night by this DIAL system are within 10 percent in the altitude 
range 30-47 km with a vertical resolution of 1 km and a horizontal resolution of 
800 km. For a case of vertical resolution of 3 km (Az=3 km, 6z»1 km) with the 
same horizontal resolution, the relative uncertainties of ozone densities will 
be reduced to 2-3 percent. The errors due to temperature dependence of the 
absorption cross section of ozone will be smaller than 1 percent since atmos- 
pheric temperature can be measured within 10 K using a wavelength of 359,8 nm, 
DIAL measurements of ozone in the lower stratosphere using this shuttle lidar 
system depend on aerosol loading and our understanding of aerosols. For the 
increased levels of stratospheric aerosols experienced after violent volcanic 
eruptions, the relative uncertainties of ozone densities will be large in the 
region below about 24 km, DIAL measurements of ozone with good accuracy will be 
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useful for Increasing our understanding of ozone in the middle and upper 
stratosphere . 

One of us (O. Ochino) is supported by NASA Grant No. NCC1-96, 
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Table I. Shuttle Lidar Parameters Assumed for DIAL Measurement of Ozone 
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Figure Captions 


I 

Fig. 1 Depiction of DIAL transmitter system. » the ith Stokes line. H 2 
is hydrogen, D 2 deuterium. 

Fig. 2 Total number of signal photoelectrons calculated by using the lidar 
system in Table 1 . 

I 

Fig. 3 Uncertainty of ozone density measured by DIAL system using for KrF laser 
fundamental and SRS of H 2 • 

Fig. 4 Uncertainty of ozone density measured by DIAL system using for KrF laser 
fundamental and SRS of D 2 • 

Fig. 5 Model extinction coefficients at a wavelength of 313.0 nm. 

Fig. 6 Uncertainty of ozone density measured by DIAL for the three aerosol 
models in Fig. 5. 

Fig. 7 Temperature variations of absorption cross sections at seven wavelengths 

20 

relative to 291 K. 

Fig. 8 Temperature model and retrievals from backscatter measurements at 359.8 

t 

j nm for two boundary conditions at 90 km: (1) Tg + 40 K, and (2) 

Tg - 60 K. 

I Fig. 9 Temperature accuracy measured at a wavelength of 359.8 nm for different 

laser output energies. Dotted line shows temperature accuracy for 
increasing quantum efficiency to 0.25 and the optical efficiency to 
0.20 for a laser output energy of 0.4 joules. 
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ATMOSPHERIC OZONE MEASUREMENTS BY EXCIMER LIDARS 

0. Uchino*, M. Maeda*, and T. Shibata* 

♦NASA Langley Research Center, Mail Stop 475 
Hampton, YA 23666, USA 

"Department of Electrical Engineering, Kyushu University, 
Hakozaki, Fukuoka 812, Japan 


Using an XeCl lidar, the first measurement of the vertical distribution of 
the stratospheric ozone was made at Fukuoka (33*N, 130*E) in July 1978 (Uchino et 
al., 1978). The method is a one-wavelength-absorption lidar with supporting 
meteorological data. Atmospheric molecular densities were obtained with daily 
meteorological radiosondes. The effect of the stratospheric aerosols on the 
measured ozone concentrations was corrected using data of scattering ratios 
observed by a ruby lidar. The method was promising for the continuous monitoring 
of the stratospheric ozone concentration (Uchino et al., 1979). 

With a more compact and reliable XeCl laser, continuous ozone monitoring has 
been made at Fukuoka since September 1979 (Uchino et al., 1980). About fifty 
data sets were obtained. Figure 1 shows a typical data of the ozone 
concentrations. At around 20-25 km the measured errors are about 10 percent. 
Above this height range, the errors increase owing to decrease of the signal. At 
an altitude of 15 km they increase up to 10-20 percent due to decrease of the 
ozone concentrations and the effect of the aerosols. 

Figure 2 shows the correlation between ozone mass mixing ratio measured by 
the XeCl lidar and temperature obtained by radiosondes at Fukuoka Meteorological 
Agency. The ozone mixing ratio has a high correlation with the temperature. 

This result is in good agreement with those obtained by ozonesondes at Tateno 
(36.1“N, 140. 1"E) (Uchino et al., 1983a). 

The effect of the stratospheric aerosol layer on the ozone concentrations 
was calculated from the scattering ratio data obtained by a frequency doubled 
Nd:YAG lidar, and estimated to be 10-20 percent in an altitude range of 15-20 km 
during 1979-1981 (Uchino et al., 1980). After the violent volcanic eruptions of 
Mt. El Chi chon in early April 1982, the concentrations of the aerosols increased 
by a factor of 100 compared with those in the background level. For these many 
and variable aerosols, we developed a high power XeCl laser pumping a dye cell of 
p-terphenyl producing two wavelengths of 340.5 nrn and 308 nm simultaneously. 

Figure 3 shows the scattering ratio profiles at 340.5 nm calculated by an 
iterative method. A parameter of A is a ratio of extinction coefficient to 
back scattering coefficient of the aerosols. The maximum scattering ratio differs 
from about 18 percent at an altitude of 25 km for three values of A (Uchino, et 
al., 1984). A value of 22.7 was used for a bimodal size distribution measured by 
an optical counter in October 1982 and assuming 75 percent sulfuric acid 
particles. The backscattering coefficient and extinction coefficient at 308 nm 
were estimated from these data using Mie theory. 


♦Visiting Senior Research Scientist from Meteorological 
Research Institute, Tsukuba, Ibaraki 305, Japan 
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Figure 4 shows the ozone concentrations measured by this method. The fluc- 
tuation of the data is large in an altitude range of 23- 26 km. It might be the 
inhomogeneous aerosols (Pelon and Megie, 1982; Browell et al., 1985) and an 
uncertainty of the optical model of the aerosols. For increased stratospheric 
aerosols, a lidar with a wavelength pair of 308 nm and 313 nm is appropriate as 
shown later. 

As the ozone concentrations decrease in the troposphere, shorter wavelengths 
are required. A high power KrF laser was developed for pumping a cell of methane 
gas. Emission of the second Stokes line (290.4 nm) of Stimulated Raman scatter- 
ing (SRS) from methane and a XeCl laser (308 nm) was used for two-wavelength 
sources (Uchino et al., 1983b). 

The ozone concentrations measured by this DIAL system are shown in Fig. 5. 
The Rayleigh component of the optical thickness was calculated from atmospheric 
molecular densities obtained by a radiosonde at Fukuoka Meteorological Agency. 

The differences of the backscattering coefficients and the Mie component of the 
optical thickness for the two wavelengths were neglected. The errors due to 
these neglects were estimated to be 1-2 percent above 4 km for a typical aerosol 
model. For more accurate estimation of errors due to the aerosols, an aerosol 
measurement at a third wavelength is available. 

In the DIAL system, the output energy of 290.4 nm was relatively low 
(2-4mJ). To obtain more efficient SRS, we further developed a high power KrF 
laser with an automatic uv preionization (Shibata et al., 1985). The output 
power was increased by a factor of two compared with the former KrF laser, and 
the energy conversion efficiencies of the Stokes lines of hydrogen gas and 
methane gas were improved as shown in Fig. 6. 

In Table 1 performances of current lidar system are shown. The maximum out- 
put energy per pulse of these SRS and excimer lasers are also shown. Figure 7 
shows the ozone profile observed on December 7, 1984. The wavelengths used to 
observe each altitude range are also shown. The measurement error are within 30 
percent except 13 - 18 km. 

As conclusions, an XeCl lidar with supporting meteorological data is useful 
for the stratospheric ozone measurement for the background aerosols. In the 
upper stratosphere, a wavelength pair of XeCl and SRS of methane or hydrogen gas 
is promising (Werner et al., 1983). 

For increased stratospheric aerosols, a lidar with a wavelength pair of 
308 nm and 313 nm is applicable although the vertical resolution is two times 
lower than that of an XeCl lidar. 

For the tropospheric ozone measurements, excimer lasers and SRS pumped by 
excimer lasers are useful. These lidar systems are comparatively simple configu- 
rations and they have possibilities of a high average power and a high repetition 
rate. 
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Table 1. Characteristics of Lidar System 
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FIGURE CAPTIONS 


Fig. 1 Ozone concentrations and their standard deviations measured by an XeCl 
lidar. 

Fig. 2 Correlation between ozone mass mixing ratio measured by an XeCI lidar and 
temperature obtained by radiosondes at Fukuoka Meteorological Agency. 

Fig. 3 Vertical profiles of scattering ratio at a wavelength of 340.5 nm for 
three values of A. 

Fig. 4 Ozone concentrations measured by two wavelengths of 308 rrni and 340.5 nm. 

Fig. 5 Ozone concentrations and their standard deviations measured by two 
wavelengths of 290.4 nm and 308 nm. 

Fig. 6 Energy conversion efficiencies of H 2 and CH 4 . 

Fig. 7 Ozone concentration measured by lidar system in Table 1. 
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